Crossover between two different Kondo couplings in side-coupled double quantum dots 
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We study the Kondo effect in side-coupled double quantum dots with particular focus on the 
crossover between two distinct singlet ground states, using the numerical renormalization group. 
The crossover occurs as the quantized energy level of the embedded dot, which is connected directly 
to the leads, is varied. In the parameter region where the embedded dot becomes almost empty 
or doubly occupied, the local moment emerging in the other dot at the side of the path for the 
current is screened via a superexchange process by the conduction electrons tunneling through the 
embedded dot. In contrast, in the other region where the embedded dot is occupied by a single 
electron, the local moment emerges also in the embedded dot, and forms a singlet bond with the 
moment in the side dot. Furthermore, we derive two different Kondo Hamiltonians for these limits 
carrying out the Schrieffer-Wolff transformation, and show that they describe the essential feature 
of the screening for each case. 

PACS numbers: 73.63.Kv, 72.15.Qm 



I. INTRODUCTION 

The Kondo effect is a prototypical many-body phe- 
nomenon that is caused by the interaction between a lo- 
calized spin and conduction electrons. Since the Kondo 
effect was observed in a quantum dot (QD) system^ 
effects of electron correlation on quantum transport have 
attracted much attention. Moreover, the recent experi- 
mental advancement enables one to examine the Kondo 
physics in a variety of systems, such as an Aharonov- 
Bohm ring with a QD and double quantum dots (DQD). 
In these systems multiple paths for electron propagation 
also affect the tunneling currents, and give rise to a Fano- 
type asymmetry in the line shapes of conductance. 

A side-coupled DQD system with a T-shape configura- 
tion, as shown in Fig. [TJ is a typical system in which the 
interplay between the Kondo effect and the interference 
effect occurs. For this type of DQD, a number of the- 
oretical studies have been carried out so far^r— It also 
becomes experimentally possible to fabricate this kind 
of geometry^ In the side-coupled DQD system shown in 
Fig- HI one of the dots (QD2), which is referred to as a side 
dot, has no direct coupling to the leads, but is coupled to 
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FIG. 1. (Color online) Side-coupled double quantum dots. 
£1(2) an d £/i(2) are the energy level and the Coulomb interac- 
tion at the QD1 and QD2, respectively, t denotes the interdot 
coupling. 



the embedded dot (QD1). Because of this unique geom- 
etry with the multiple paths, intriguing phenomena can 
occur in this system. For example, a two-stage Kondo 
effect can occur for the small interdot coupling t, where 
each of the dots is occupied by a single electron.— ^ 12 ' 14 ' 15 
In this case, the local moment in the QD1 is screened first 
by the conduction electrons at higher temperature, and 
then the moment in the QD2 is screened at lower tem- 
perature to form a singlet ground state. Most of the pre- 
ceding studies of a side-coupled DQD system focused on 
the two-stage Kondo effect taking place in this situation. 
In contrast, for large t, the two adjacent local moments 
screen each other to form a molecular-type singlet. 

The gate voltage that is applied to the dots can further 
change the charge and spin distributions in the DQD, 
and evolve the system toward the mixed- valence regime. 
Specifically, as the energy level E\ of the QD1 is varied, 
another typical singlet state appears in the parameter 
region where the QD1 becomes almost empty or dou- 
bly occupied. It is a singlet bond between the local mo- 
ment at QD2 and the conduction electrons, and is formed 
by a superexchange mechanism. Some numerical indica- 
tions that this type of singlet state is formed were seen in 
data of previous works.^ Maruyama et al. obtained an 
asymmetric conductance peak of the Fano shape for fi- 
nite Ei in the Kondo regime^ However, their study is 
focused mainly on the transport properties. Zitko et 
al£ examined a similar situation, and showed that the 
Kondo temperature becomes small in this case. The pre- 
cise features of the singlet state, however, were not fully 
examined. Some groups have studied the gate-voltage 
dependence of the conductance where the energy levels 
of the two dots are moved simultaneously.— However, 
it has still not been clarified in detail how the singlet 
ground state evolves across the crossover region between 
the singlet state due to the two-stage Kondo effect and 
the one formed by the superexchange mechanism between 
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the QD2 and the conduction electrons. 

In this paper, we re-examine a side-coupled DQD sys- 
tem shown in Fig. [1] and study how the Kondo singlet 
bond is deformed as the energy level e\ in the QD1 is 
changed. We find that as E\ moves away from the Fermi 
energy, the electrons at the QD1 cannot contribute to the 
screening of the local moment at the QD2, and the con- 
duction electrons tunneling into the QD2 virtually via e\ 
screen the local moment. This electron tunneling process 
is similar to the superexchange mechanism seen in tran- 
sition metal oxides such as MnO and CuOi 17 ' 18 In the 
present case, the Kondo singlet bond is formed between 
the QD2 and the leads, mediated by the quantized level 
e\ at the QD1. This mechanism is quite different from 
the Kondo screening in the case of e\ ~ where a sin- 
glet bond between the QD1 and QD2 plays a dominant 
role.— ~— Therefore, the gate voltage applied to the QD1 
deforms the Kondo cloud, and it can be probed through 
the variation in the phase shift of the DQD. In order 
to clarify these features, we calculate the phase shift us- 
ing the numerical rcnormalization group. Furthermore, 
we calculate the spin susceptibility to obtain the Kondo 
temperature, which shows good agreement with the one 
obtained from an effective Kondo Hamiltonian we have 
derived in this work. In the presence of the Coulomb 
interaction U\ at the QD1, we also find that the two- 
stage Kondo screening changes to a single-stage process 
as £\ moves away from the electron- hole symmetric point 
ei ~ -E/i/2- 

This paper is organized as follows. In Sec. [TTl we give 
the Hamiltonian of our system. In Sec. IIII1 we derive 
the effective Hamiltonian using the perturbation theory 
in the tunneling matrix elements, which is identical to 
that derived from the Schrieffer- Wolff transformation.— 
In Sec. IIV1 we show the numerical results, and discuss 
how the Kondo singlet state evolves as e\ varies. A sum- 
mary and discussions are given in Sec. [V] 

II. MODEL 

The Hamiltonian of a side-coupled DQD system shown 
in Fig. Q] reads 

H = Hqdi + Hqd2 + Hint + ^ (H v + Ht, v ), (1) 
where 

(7 

H in t = t ^ (d\a d 2a + H - C -) , H » = Y £ ^ka^.far , 
a k.cr 

H ^ = E ( c U d ia + H - c -) , » = L,R. (2) 

Hqdi describes the QD1 for i = 1, and the QD2 for 
i = 2, £, the energy level, Ui the Coulomb interaction, 



and nj jCT = d ia d ir7 . H mt denotes the interdot coupling 
with the hopping matrix element t. H L / R describes the 
normal lead of the left/right side. Vz,m is the tunneling 
matrix element between the QD1 and the left/right lead. 
We assume that T L / R (e) = nVl, R J2k — £ k)/N is a 
constant independent of the energy e, where N is the 
number of the states in each lead. 

For convenience of the following discussions, we apply 
a unitary transformation to the leads, using the inversion 
symmetry, 

_ VhCL.ka + Vrcrm _ Vlclm - Vrcrm 

Ska — i a ka -rr 5 

V s Vs 

V s = ^Jv L 2 + V R 2 . (3) 

Then, the Hamiltonian for the leads, Yli>e{L R}(Hv + 
Ht,v), can be rewritten as 

{H v +H T>v ) = H s + H a + H TiS , (4) 

v£{L,R} 

with 

^, S = E^( S L^+H.c). (6) 

Note that the operator Sk a couples to the QD1, while 
the operator ata is decoupled. Therefore, we can map 
the original model given by Eq. ([T]) to the two-impurity 
Anderson model (TIAM), 

%TIAM = Hqdi + Hqd2 + Hi n t + H s + Ht.s- (7) 

Using this model, we will discuss how the singlet state 
due to the Kondo effect changes as the energy level £\ at 
the QD1 varies. 

III. KONDO HAMILTONIANS 

In this section, we derive the effective Hamiltonians 
in order to study the Kondo behavior in the two oppo- 
site cases, £\ ~ and large e\. For this purpose, we 
use the perturbation theory in the tunneling matrix el- 
ements t and Vs, which is equivalent to the Schricffcr- 
Wolff transformation^ We assume that the Coulomb 
interaction at the QD1 is zero, U\ = 0, in this section 
in order to focus on the effects of e\. We also consider 
the effects of U\ in Sec. IIVB and the Appendix. 

Figure [2] shows a schematic energy diagram of a side- 
coupled DQD for £2 + ^2/2 — 0. Let us first consider 
the case, shown in Fig. [Ha), where the energy level at 
the QD1 is located around the Fermi energy e\ ~ 0. 
When the Coulomb interaction U2 is much larger than 
the interdot coupling t, the local spin moment arises at 
the QD2. Carrying out the perturbation expansion in 
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FIG. 2. (Color online) Schematic energy diagram of a side- 
coupled DQD. (a) £1 ~ 0, and (b) ei > such that (ni jCT ) ~ 
is satisfied. Ep denotes the Fermi energy of the leads. 



the tunneling matrix element i, we obtain the Kondo 
Hamiltonian of the form 



U^k - J12 ^2 d\ a -^-d l0 



where 



QDl 



H, + H 7 



(8) 



2tr 
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n=0.e2 = -U2/2 8t 

~ >U a 



£i {s2 + U 2 )-ei 

2 



(9) 



T aa , — (t* c , , T% a , , t* c , ) is the vector representation of 
the Pauli matrix, and S2 is the spin operator of the local 
spin at the QD2. In the electron- hole symmetric case, 
E\ = and e 2 = —U2/2, the exchange coupling is given 
by J12 = 8t 2 /U2- This type of exchange coupling ap- 
pears in the usual Kondo model, and is referred to as 
the "interdot exchange (IE)" coupling in the following, 
in order to distinguish it from the superexchange cou- 
pling discussed later. In the case where J12 is larger 
than r (= Tl + ^r), the local moments at the QDl and 
QD2 form a molecular-type singlet bond. In the opposite 
case J12 < r where the QDl is coupled more strongly to 
the leads, the spectral weight of the electron at the QDl 
is broadened, and the local moment S 2 at the QD2 is 
screened by the electrons with this broadened density of 
states. In both of these cases, the singlet state is formed 
mainly between the moments in the QDl and QD2. Sim- 
ilar considerations also make sense for finite Uir& just 
by replacing T with the Kondo temperature T^ D1 for 
the QDl. Namely, the two-stage Kondo effect occurs 
for J \2 < T® D1 , whereas a singlet bond becomes the 

molecular-type one for J\ 2 > T^ 1 . 

Next, we consider the situation shown in Fig. EJb), 
where the energy level at the QDl is away from the Fermi 
energy. The Kondo effect in this situation is a main focus 
of this paper. In this case, the QDl is almost empty (dou- 
bly occupied) for e\ > (ei < 0), and the spin degree of 
freedom disappears at the QDl. Thus, the electrons at 
the QDl can not contribute to the screening of the local 
moment at the QD2, and E\ works as a potential barrier 
that disturbs charge transfer between the QD2 and the 



leads. Maruyama et al^ and Zitko et al£ also consid- 
ered the |ei| — > 00 limit. However, the precise features 
of the screening process due to the superexchange mech- 
anism were not examined in detail. A similar situation 
to Fig. [UJb) also arises in transition metal oxides such as 
MnO and CuO, where the superexchange interaction de- 
scribes the coupling between the local spins in magnetic 
ions mediated by nonmagnetic oxygen anions i ' 18 

In order to clarify the Kondo screening in the situa- 
tion shown in Fig. d{b), we derive the effective Hamil- 
tonian from the perturbation expansion with respect to 
the tunneling elements t and V s to the fourth order (see 
the Appendix). The result can be expressed in the form 
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where 
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The coupling constant Jse depends on the energy level 
Ei, which is caused by a virtual process with a single 
electron passing through the QDl. This term appears 
as the Vgt 2 -type contribution in the fourth order per- 
turbation expansion with respect to tunneling elements. 
The screening of the local moment S2 is achieved for 
large £1 by the conduction electrons tunneling virtually 
through the QDl. This screening mechanism is essen- 
tially the same as the one due to the superexchange in- 
teraction mentioned above, and the singlet bond becomes 
long compared to that in the case of E\ ~ 0. Therefore, 
Jse is referred to as the "superexchange (SE)" Kondo 
coupling in the following. Note that a similar screening 
occurs also for negative £i(< 0), although Fig. [^b) de- 
scribes only the situation for positive £i(> 0). For nega- 
tive large £1, the QDl is almost doubly occupied and the 
effective Hamiltonian takes the same form ?4 b) in Eq. 

(USD. 



IV. NRG RESULTS 

A. Crossover between the IE and SE Kondo 
screenings 

To confirm the above discussions more precisely, we 
calculate the phase shift, the average number of electrons 
in each of the dots, and the spin susceptibility for the 
two-impurity Anderson model "Htiam using the numeri- 
cal renormalization group (NRG)?22, We first show the nu- 
merical results of the phase shift due to the DQD, which 
is helpful to clarify the formation of the singlet state be- 
cause the phase shift reflects an electron scattering at the 
DQD. From the phase shift <p, we can also deduce the to- 
tal number of electrons Noqd = X)<r( n i-°' + n 2.a) in the 
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DQD, using the Fricdcl sum ruL 



N 



DQD 
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Figure [3] shows Ndqd as a function of Ei/t, and the num- 
ber of electrons in each of the dots, (^1(2)) = Tl a ( n l{2),a) ■ 
Let us first look at a region around e± = 0. In this re- 
gion, namely Noqd — 2, both of the dots are nearly 
half-filled ((ni(2)) — 1), and the phase shift takes the 
value of if ~ 7r from Eq. (|12jl . The singlet state is formed 
dominantly inside the DQD, and thereby the conduction 
electrons at the leads are not scattered by the local spin 
at the dots. As e\ increases, Ndqd shows a sharp drop 
around £\jt = 0.8 and approaches Njjqd — 1- From Eq. 
(|12p. we see that the phase shift if also changes from 7r 
to ir/2. In the region around e\jt = 2.0, n\ goes to zero 
while n 2 almost remains unchanged. It indicates that 
the local spin appears only at the QD2 for large £i, and 
the SE Kondo coupling can be described by the Kondo 
Hamiltonian given in Eq. (|10[) . Therefore, the kink be- 
havior of the phase shift with the height 7r/2, seen in Fig. 
131 signifies the crossover between the IE Kondo screening 

fa") 

described by the Hamiltonian with J\ 2 and the SE 
one described by H$ with Jse- 

In order to discuss the change of the low-energy states 
for the energy level e± in more detail, it is helpful to use 
the fixed-point Hamiltonian in terms of the renormalized 
parameter; 



(13) 



,2:s 



H$ = e 2 n 2 + emx + t ^ (d\ a d 2a + H.c. 
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FIG. 3. (Color online) Total number of electrons in the DQD, 
Ndqd, and each number of the electrons at the dots. We set 
Ui = 0, U 2 /t = 6, e 2 = -t/2/2, and T L /t = F R /t = 0.1. The 
inserted figures illustrate a formation of a dominant singlet 
bond which is described by the (green) dashed line. 




FIG. 4. (Color online) (a) Renormalized parameters £2, t and 
(b) conductance at zero temperature as a function of ei/t. 
The parameters of the system are the same as in Fig. [3] Inset 
of (a): Phase boundary between singlet and doublet ground 
states for the isolated DQD system (Tl/r = 0), where we set 
Ui =0 and £2 = -1/2/2. 



where 



e 2 = Z{e 2 + S 2 (0)), t = Vzt, 



Z= 1- 



3S 2 (e) 



ds 



(14) 



S 2 (e) is the self energy due to the Coulomb interaction 
C/ 2 . Using this fixed-point Hamiltonian, a unified analysis 
for the £1 dependence becomes possible. The crossover 
between the two opposite limits, at Si = and E\ — > 00, 
can be described as a continuous change of the parameter 
values of the fixed-point Hamiltonian. We can calculate 
the renormalized parameters e 2 and t using the NRG^ 
The results are shown in Fig.|U[a). It is noteworthy that 
both £2 and t show a sharp decrease around £\jt ~ 0.8, 
which indicates the crossover between the two different 
singlet bonds, namely the one due to the IE coupling 
J12 and the other due to the SE Kondo coupling Jse- 
Furthermore, from these parameters we can deduce the 
phase shift if of the DQD and the conductance G at zero 
temperature £ 



if = — + tan 



2e 2 



t - £l£\ 



e 2 r 



2e 2 



G = -sin> = -<;i 



t 1 - £ 2 £i 

e 2 r 



(15) 



(16) 



where T = + Dr. Note that this expression for G 
is exact at zero temperature for the symmetric coupling 
= r#, and can be obtained, for instance, by using 
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the Meir-Wingreen formula^ for the Hamiltonian in Eq. 
(P) . The phase shift <p and the conductance G can be de- 
duced from the exact NRG results for the renormalizcd 
parameters. Figure H^b) shows the result of the conduc- 
tance as a function of E\ for £2 = —U2/2. We see that 
the conductance shows an upturn around £i/t ~ 0.8, and 
at this value the crossover between the ground state due 
to the IE coupling and that due to the SE Kondo cou- 
pling occurs. The behavior of the conductance in Fig. 
@Jb) can also be explained in terms of the Fano-Kondo 
effect. This is because the energy level £1 at the QD1 
varies the asymmetric parameter q for a Fano line shape, 
as discussed by Maruyama et al~ and ZitkoJ^ Indeed, 
the conductance decreases at e\/t ~ 0, where q ~ 0, due 
to the destructive interference effect while the conduc- 
tance approaches 2e 2 /h in the limit of £\/t — > 00 where 
q — > 00. 

The nature of the crossover can also be related to a 
level crossing taking place in a molecule limit Tl/r = 0, 
where the QD1 is decoupled from the lead. In this 
limit the isolated DQD is described by a Hamiltonian 
HQDi+HQ02+Hi nt , and the ground state of the molecule 
becomes a singlet or doublet, depending on the value of 
\ei\/t and U$lt. As shown in the inset of Fig. Q|a) for 
U\ = 0, the ground state is a spin singlet if either \ei\/t 
or L^/i is small. In the opposite case, a spin doublet be- 
comes the ground state. Note that in the doublet region, 
nearly one electron occupies the QD2, whereas the QD1 
is almost empty or doubly occupied. Thus, the local mo- 
ment emerges mainly at the QD2. We see in the phase 
diagram in Fig. 0Ja) that the transition takes place in 
this molecule limit at |£i|/£ — 0.8 for t/2/t = 6, and it 
agrees well with the position where £2 and t show a sharp 
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decrease. For finite V 



L/R: 



the conduction electrons can 



tunnel from the lead to the QD2 via the QD1. However, 
the electrons at the QD1 cannot contribute to the screen- 
ing of the moment at the QD2, because the QD1 is almost 
empty or doubly occupied, and has no local spin moment. 
Then, the local spin at the QD2 is screened by the con- 
duction electrons from the leads over the QD1, which is 
the SE Kondo screening discussed above. Therefore, us- 
ing the phase diagram in Fig. E{a), we can estimate the 
value |£i|/t and Uz/t, at which a crossover between two 
distinct singlet states occurs. 

In order to estimate the Kondo temperature due to the 
SE process, we calculate the contribution of the QD2 to 
the impurity susceptibility, defined by 
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FIG. 5. (Color online) Plots of k B T X 2/(gfJ.B) 2 vs k B T/D 
for several values of ei/t. We set L^/i = 6, £2 = — U2/2, 
r L /t = F R /t = 0.1, and t/D = 10~ 3 . The arrows with dashed 
lines indicate the Kondo temperature Tk(Jse) obtained from 
Eq. Q5J| for ei/t = 1.2, 1.6, and 2.0 (see Table QJ. 
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Tk(Jse) = D^J^ cxp(-l/V SE ), (18) 
where D is the half-bandwidth of the leads and p = 1/2D. 

TABLE I. The Kondo temperatures T K (J SE ) from Eq. (JTHJl. 



ei/t 


1.2 


1.6 


2.0 


T k (Jsb)/D 
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1.45 x 10~ 14 
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The Kondo temperature Tk{Jse), which is obtained 
from Eq. (fT5|) , is listed in Table Q] for several values of 
E\jt. Furthermore, these values of Tr-(Jse) are indicated 
by the arrows with dashed lines in Fig. [5j We see that 
T X 2 obtained from the NRG for e x jt = 1.2, 1.6, and 2.0 
decreases rapidly showing a clear crossover to the Kondo 
regime around the temperature indicated by the arrows. 
This agreement demonstrates that the Kondo screening 
for a large value of £1 is mainly owing to the SE one 
described by the Kondo Hamiltonian . 



B. Coulomb interaction Ui at QD1 



X2 



(ffMfi) 2 



{(Si) - (S 2 Z ) ) 



(17) 



where (S 2 ) ((S%)o) is the z component of the total spin 
of the system with (without) the QD2. Figure [5] shows 
the results of \ 2 for several values of e\. We can es- 
timate the Kondo temperature from the slope of these 
plots In particular, for a large value of £1, we can com- 
pare the Kondo temperature estimated from T\2 with 
that obtained from the Kondo Hamiltonian Hv with a 



So far, we have assumed that the QD1 is noninteract- 
ing, U\ = 0. In this subsection, we discuss the effects of 
the Coulomb interaction C/i at the QD1 on the energy 
scale of the Kondo screening. 



1. Finite Ui for Ei > 

We first introduce the Coulomb interaction Ui for pos- 
itive £1 > 0. Figure O shows the impurity susceptibility 
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FIG. 6. (Color online) Plots of k B Tx2/{g^ B f vs k B T/D for 
several values of ei/t and Ui/t, where we set Uz/t = 6, £2 = 
-U2/2, T L /t = r R /t = 0.1, and t/D = 1(T 3 . 



at the QD2, \2, for several values of E\ and U±. The other 
parameters are the same as those of Fig. [5J Since we set 
£1 in a way such that there is almost no local spin mo- 
ment at the QD1, the SE Kondo screening occurs in these 
examples. In Fig. [51 we see that the crossover tempera- 
ture, corresponding to the Kondo energy scale, decreases 
as XJ\ increases. This is because the QD1 becomes almost 
empty for large positive £\ and large Coulomb repulsion 
\J\, and the virtual electron tunneling from the leads to 
QD2 is suppressed significantly. Specifically, in the case 
where E\ is away from the Fermi energy {s\ > t,T), the 
SE Kondo coupling Jse between the QD2 and the leads 
given in Eq. (jll[) can be expressed in a more general 
form, by taking into account the effect of the Coulomb 
repulsion as 



Jse 



Vst 

£1 



1 



g?C/i - (e 2 - 2 £l )( £l - e 2 ) 2 | 
e 2 (£2 - 2ei - J7i)(£i - e 2 ) 2 £2 + ^2 



(19) 



This Jse monotonically decreases with increasing Ux, so 
that the Kondo temperature Tk(Jse) defined by Eq. (fTS)) 
also becomes small as U\ increases. We can see the cor- 
responding shift of the crossover temperature in Fig. [6] 
for the two-impurity Anderson model 'Htiam- 



Crossover between the two-stage and sing 
Kondo screenings 



-stage 



Next we consider another case: the crossover between 
the two-stage Kondo screening and a single-stage Kondo 
screening. It has been discussed previously that the two- 
stage Kondo effect can occur in the case where each of the 
two dots has a local moment and the Kondo temperature 
T® D1 for the QD1 is larger than the exchange coupling 
J12 between the dots£i2 It takes place typically near the 
electron-hole symmetric point, where £1(2) + L/i(2)/2 — 0. 
Therefore, as the energy level of the QD1 moves away 
from the electron-hole symmetric point E\ ~ —Ui/2, the 
SE Kondo screening can arise because the local spin mo- 
ment at the QD1 disappears. Furthermore, the two-stage 
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2 , (b) k B Tx2/(g^B) 2 , and (c) entropy 



FIG. 7. 

SuQu/ks as a function of k B T/D for several values of ei/T 
We set t/T = 0.3, Ui/T = U 3 /T = 6, e 2 /U 2 = -0.5, and 
T/D = 10~ 3 . 



Kondo screening near the symmetric point changes to the 
single-stage Kondo screening. This was also discussed 
partly by Zitko et al. They showed that the Kondo 
temperature rapidly drops as E\ moves away from the 
electron-hole symmetric points However, how the sin- 
glet ground state evolves in the crossover region has not 
been clarified in detail. In order to confirm the precise 
features of the Kondo screening, we calculate the suscep- 
tibility for the two dots as well as that for the QD2, 

XDQD = { -^f - (S^cad) • (20) 

Here, (S^) is the z component of the total spin of the 
whole system including the DQD, and (S^)\ e ad is the 
same quantity without the DQD. Furthermore, we also 
calculate the entropy of the DQD defined by 



(21) 



Here, E = TT[He- H ^ kBT ^]/Ti[e- H ^ kBT ^\ and F = 
— fcsrinTr[e^ H ^ fefjT '] are the internal energy and free 
energy of the whole system consisting of the DQD and 
the leads, while i?i ea d an d -Fi ea d are those for the uncon- 
nected leads. 
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FIG. 8. (Color online) The results presented in Fig. [7] are 
replotted over a wide temperature range. These results con- 
firm that the moment is eventually screened in the cases of 
£i/r = 1 and 2 although the screening temperature becomes 
exponentially small. Therefore, to observe the screening pro- 
cess due to the SE Kondo coupling Jse defined by Eq. (|19[) 
in a realistic temperature, e\ should be not so large. 

Figure [TJa) shows the spin susceptibility xdqd for sev- 
eral values of £i/r. For £i/T = —3, which corresponds to 
the electron-hole symmetric point E\ + U\/2 = 0, Txdqd 
shows a peak around T/D ~ 10~ 3 . This peak indicates 
that each of the dots is occupied by a single electron, and 
the local moment is well developed. In this case, a two- 
stage screening occurs as temperature decreases: the first 
stage can be seen at T/D ~ 10~ 4 , and the second one 
at T/D ~ 1CT 7 . These two energy scales correspond to 
the Kondo temperature for the first stage and that 
for the second stage T|- nd , respectively. If the Coulomb 
interaction at the dots is much larger than the tunneling 
constants T and t, then the peak of T^dqd approaches 
0.5 in units of {g^s) 2 and the structure that emerges at 
Tjf" becomes clear. For £i/T = —2 and — f, which are 
still negative but closer to the Fermi energy, the two- 
stage Kondo effect can be seen more clearly. The Kondo 
temperature T|- nd for the second stage becomes lower: 
T% nd /D ~ 1CT 8 for £i/T = -2 and T$? d /D ~ ICT 11 for 
£i/r = —1. In Fig. [3[c), we see that the two-stage be- 
havior can be observed more sharply in the temperature 
dependence of the entropy than that of the susceptibility. 

As Si crosses the Fermi energy and takes a positive 
value, the peak of Txdqd seen in Fig. [TJa) at T/D ~ 



10 -3 is suppressed. This indicates that the local mo- 
ment at the QD1 disappears as the QD1 becomes almost 
empty. Then, the local moment at the QD2 is screened 
by the SE process via a single stage as the temperature 
approaches zero. We see in Fig. [TJb) that T\2 is al- 
most constant (0.25) in a wide temperature range below 
T/D~I0~ 3 . This indicates that the local spin moment 
at the QD2 remains almost free in this temperature re- 
gion. Furthermore, the entropy for the DQD shown in 
Fig. He) is locked at the value of ln2 below T/D ~ 10~ 4 , 
which is caused by the unscreened spin moment at the 
QD2. 

This free moment must be screened eventually at low 
temperatures by the SE Kondo coupling Jse although the 
Coulomb interaction U\ makes the Kondo temperature 
very small. In order to see the low temperature region, 
Xdqd, X2, and Sdqd are shown in a wide temperature 
range in Fig. [5] We see that Txdqd for £i/r = 1 and 
£i/r = 2 shows the decrease at T/D ~ 10" 45 and T/D ~ 
10 -90 , respectively. These results confirm clearly that 
the local moment at the QD2 is really screened although 
the energy scale for the SE Kondo screening becomes 
small for large U\ and £\. 



V. SUMMARY AND DISCUSSION 

We have studied the Kondo effect in a side-coupled 
DQD system with focus on how the Kondo singlet state 
changes by varying the energy level at the embedded dot 
(QDI ). We have found that when the side dot (QD2) is in 
the Kondo regime, two distinct singlet states appear; one 
is due to the IE coupling between the QDI and the QD2, 
and the other is caused by the SE Kondo coupling be- 
tween the QD2 and the leads via the QDI. In this sense, 
the latter is a different type of the singlet state from 
the former which has been studied in the side-coupled 
DQD systems so far. In order to clarify the screening 
process, we have obtained the effective Kondo Hamilto- 
nians using the perturbation expansion with respect to 
the tunneling matrix elements. From these Kondo Hamil- 
tonians, we have shown that in the case where the QDI 
is almost empty and doubly occupied the screening is 
caused by a superexchange mechanism, and the singlet 
bond becomes long. Moreover, we have calculated the 
phase shift and the conductance using the NRG method, 
and have obtained the relation between the phase shift 
and the conductance. We have found that the conduc- 
tance is enhanced at the crossover region between the 
singlet ground state due to the IE coupling and that due 
to the SE Kondo coupling. We have also calculated the 
local spin susceptibility and have estimated the Kondo 
temperature, which shows good agreement with that ob- 
tained from the effective Kondo Hamiltonian. Further- 
more, we have demonstrated precisely how the two-stage 
Kondo screening changes to a single-stage process as the 
energy level at the QDI moves away from the electron- 
hole symmetric point. 
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In closing, we would like to make some comments on 
the SE Kondo screening. The scenario of the SE Kondo 
screening is not limited to our side-coupled DQD sys- 
tem, but is more generic for the Anderson model where 
the impurity spin and the conduction electrons are con- 
nected via the discrete energy level. In usual cases of the 
Kondo problem, a magnetic impurity is coupled directly 
to the conduction electrons, for instance, a bulk system 
with a magnetic impurity. On the other hand, in the 
two-impurity Anderson model we have considered, a local 
spin moment is coupled indirectly to the conduction elec- 
trons via a discrete energy level. The resulting screening 
process shows a unique feature, which has clearly been 
demonstrated in this paper with the help of the effec- 
tive Kondo Hamiltonians and the NRG method. The SE 
Kondo screening discussed in this paper also appears in 
a side-coupled DQD system coupled to normal and su- 
perconducting leads, which we have previously studied^ 
In this case, the superconducting proximity to the em- 
bedded dot quenches the local spin moment because this 
proximity tends to make a singlet consisting of a linear 
combination of the empty and doubly occupied states. 
Thus, the superconducting proximity to the embedded 
dot plays the role of a potential barrier between the side 
dot and the normal lead, which can cause the SE Kondo 
screening. 

Recently, the side-coupled DQD system has been fab- 
ricated in experiments.— We thus expect that in the near 
future, it may become possible to observe the SE Kondo 
screening discussed in this paper, providing further in- 
teresting examples of correlation effects in the context of 
electron transport in nanoscale systems. 



We thus choose the unperturbed Hamiltonian to be 

H = H QD1 +H QD2 + H S . (22) 

The unperturbed ground state is described by d\ a \F), 
where \F) is the Fermi sea of the conduction band. The 
tunneling terms Ht, s and Hi nt are taken to be as the 
perturbation Hamiltonian, as illustrated in Fig. [9l 



QD2 QD1 Lead 




FIG. 9. (Color online) Schematic energy diagram of a side- 
coupled DQD for si > 0. This situation is the same as that 
for Fig. [U^b). In this figure, we describe the virtual tunneling 
events by the (blue) dashed arrows, which are labeled by a, b, c 
and d. The bold arrow at the QD2 indicates the local spin. 

The supcrcxchangc mechanism via the energy level ei 
at the QD1 is described by the fourth order perturbation 
with respect to the tunneling matrix elements t and V s . 
In the fourth order perturbation, there are six processes. 
Using the label for virtual electron tunnclings shown in 
Fig. [5J these six processes are given by 
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APPENDIX: DERIVATION OF THE EFFECTIVE 
KONDO HAMILTONIAN 



We outline the derivation of the effective Hamiltonian 
for the SE Kondo coupling given in Eq. (|10[) by using the 
perturbation theory in the tunneling matrix elements.— 
The unperturbed ground state is chosen to be the one 
for t = V s = and -e 2 , e 2 + U 2l £i > E F (= 0). In this 
situation, the ground state of the two-impurity Anderson 
model defined in Eq. ([7]) is described by the singly occu- 
pied state at the QD2 and the empty one at the QD1. 



(i) a -> b -> c -> d 

(ii) a ->• c^r b ->• d 

(iii) a —t c — >• d — > b 



(iv) c^rf-)o^() 

(v) c -> a -> b -> d 

(vi) c->a-)rf-)(). 



The virtual electron tunnclings described by the labels 
a and d are caused by H lnt , whereas those described by 
labels b and c are caused by H T s . Therefore, the fourth 
order perturbation takes the form 



1 



1 



1 



-Ha, 



E — Hq E — Hq E — Mq 



(23) 



where Hj for j G {1,2,3,4} depends on the virtual pro- 
cesses, namely Hj = H mt for the virtual electron tun- 
nelings labeled a and d, whereas ifj = Ht, s for those 
labeled b and c. As a representation of Eq. (f2U| in the 
subspacc of d\ a \F), we obtain the effective Hamiltonian 
in the form 



= ^ A E E -L-va + j x E E V-^fsv* ■ 

(24) 



fc.fc' a 



where A runs over (i) to (vi). The couplings W x and J* 
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are given by 



2e 2 (e 2 -£i) 2: 
{V s tf 



2(e 2 -2e 1 )(e 2 -e 1 ) 2 



J«=4W (i) , (25) 
, = (26) 



(Vstf 



2ei(e 2 - 2ei)(e 2 - ei) 



j( iii ) = -4VK (iii) , 



2(ei) 2 (e 2 + C/ 2 ) 

^(v) =W -(iU) j j(v) = J (iii) 5 



2(e 2 -2e 1 )(e 1 ) 2 ' 



(27) 

, j( iv ) = -41T< iv ), (28) 
(29) 

j( vi ) = _4W(vi). (30) 



Summing up these six processes, we obtain the effective 
Hamiltonian for the SE Kondo coupling, corresponding 



to ?4 b) given in Eq. ([TO]). 

For the finite Coulomb interaction U\ at the QD1, the 
processes of (ii), (iii), (v), and (vi) depend on U\ since 
the QD1 is doubly occupied in the intermediate state. 
Thus for finite U\, W x takes the form 



00- 



w 



2{e 2 -2e 1 -U 1 ){e 2 -e l f 
{Vstf 



(31) 
(32) 



2ei(e 2 -2£i- E/i)(e 2 -£i) : 
W^=W {m \ (33) 

wW- {Vst)2 



2(e 2 -2e 1 -C/ 1 )(e 1 ) i 



(34) 



Correspondingly, J A for A = (ii), (iii), (v), and (vi) is 
given by J x = —4W X , and thus we can obtain the ex- 
change coupling Jse given in Eq. ([T9|) . 
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